Estimations for the Higgs boson production with QCD and EW 
corrections in exclusive events at the LHC 



O 

(N 



Oh 
C/3 



c3 



M. B. Gay Ducati 
High Energy Physics Phenomenology Group, GFPAE, 
Instituto de Fisica, Universidade Federal do Rio Grande do Sul, 
Caixa Postal 15051, CEP 91501-970, Porto Alegre, RS, Brazil 



G. G. Silveira 

High Energy Physics Phenomenology Group, GFPAE, 
Instituto de Fisica, Universidade Federal do Rio Grande do Sul, 
Caixa Postal 15051, CEP 91501-970, Porto Alegre, RS, Brazil and 
Q_i. Center for Particle Physics and Phenomenology (CPS), 



Universite catholique de Louvain, B-13^8 Louvain-la-Neuve, Belgium 

Abstract 



The Higgs boson production is investigated in proton-proton collisions at next-to-leading-order 

00 ' 

■ accuracy in central exclusive diffractive processes at the LHC. The production process by the dou- 

ble Pomeron exchange is analyzed in the diffractive factorization through the Ingelman-Schlein 
\ approach, taking into account the parton content of the Pomeron by the diffractive partonic dis- 

tribution function provided by the HI Collaboration. Hence, we estimate the production cross 
section of the Higgs boson as well as its rapidity distribution for distinct energies of the LHC. Also, 
we include the gap survival probability in our calculation, which is studied in recent works and 



expected to lie in the range between 1% and 5% for the energy regime of 14 TeV. As a result, we 
found a production cross section of about 0.3-0.8 (1.2-3.7) fb at 7 (14) TeV, being of the same 
order as predicted by the two-photon and the Balitsky-Fadin-Kuraev-Lipatov Pomeron mecha- 
nisms. Therefore, assuming the selection rules of spin-parity properties, the exclusive production 
is a promising channel for the Higgs boson detection in the LHC. 

PACS numbers: 12.15.Lk, 12.38.BX, 13.85.Dz, 14.80.Bn 
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I. INTRODUCTION 



The inclusive processes provide the highest cross sections for particle production in hadron 
colliders, which is expected to be the discovery channel of the Higgs boson at the LHC 
Particularly, the highest contribution for this production mechanism is the gluon fusion 
vertex pp — > (gg) — >■ H + X [2], predicting a total cross section of about 15 (50) pb at y/s = 
7 (14) TeV for a Higgs boson mass of Mh = 120 GeV for different methods for the next-to- 
next-to-leading-order (NNLO) perturbative calculation [3]. Moreover, there are processes 
where the final state is the same as the decay products of the Higgs boson, that attenuate 
the production signal, e.g., the processes gg ^ bb 0], especially for a Higgs boson mass 
Mh < 135 GeV, and the gg — > 77 through the box diagram [5]. So, these events have to 
be suppressed in order to enable the Higgs boson detection, and higher-order perturbative 
calculations are highly necessary. 

On the other hand, the central exclusive diffractive (CED) process pp — > p + [LRG] + 
H + \LRG] + p [6fl is a rich tool to investigate particle production in a cleaner environment 

Q fl fl 

and was already studied for the production of x mesons [7j], dijets [8|, diphotons [9|, and the 



Higgs boson [lOj. The main signature of this process is the large rapidity gaps present in 
the final state, i.e., a region in the rapidity distribution with no hadronic activity, since the 
interaction by the exchange of Pomerons occurs with no change of the quantum numbers 



of the particles in collision, with the protons remaining intact in the final state [ill ]. An 
advantage of this process is the possibility to suppress high background signals by means 
of the J z selection rule [l2], that increases the signal-to-background (S/B) ratio for a 
Higgs boson detection at the LHC. For instance, the exclusive gg — > bb process has a cross 
section of the same order as the exclusive Higgs boson production at the LHC, and then 



the background process may overcome the production signal 



131 ] . Thus, some detectors 



are going to be set up in the LHC experiments to detect the rapidity gaps, increasing the 
possibility of observation of such processes [jjj]. 

In addition, secondary interactions occur between the protons during the exclusive reac- 
tion that causes a contamination of the final state by other particles, reducing the signal of 
the CED process. In order to weigh the fraction of processes where the secondary interac- 
tions (underlying events) do not fill the rapidity gaps, one has to compute the rapidity gap 
survival probability (GSP), which is used to obtain the cross section that will be observed 



experimentally. The studies performed by experimental groups have shown that the use of 
an overall suppression factor is favored by the phenomenological analysis with the data of 
diffractive dijet production from the Tevatron experiments 15f . There are different groups 
computing the GSP for the CED Higgs boson production 
1% and 5% is expected for the LHC kinematical regime [l8|, Il9| . 

There are two different approaches to compute the CED Higgs boson production that 



JA 1,113, and a probability between 



he first one consists in the calculation based 



10[ to compute the scattering amplitude of 



depend on the dynamics of the hard Pomeron. 
on nonperturbative [201] or perturbative QCD 
two-gluon exchange in the t-channel. In the latter, the diffractive interaction is set in a 
way that the gluon fusion vertex is preceded by a soft gluon exchange in order to neutralize 
the color flow into the bosonic loop, introducing the essence of the Balitsky-Fadin-Kuraev- 
Lipatov Pomeron: an exchange of two hard gluons in the t channel. Additionally, some 
phenomenological aspects are also included to consider important features of this process, 
like off-diagonal unintegrated gluon density, and Sudakov form factors at leading logarithm 
approximation, that have been recently reexamined 21]. 

The second possibility, in which we are mainly interested in this work, was proposed in 
an analysis of the jet structure in the LEP data [22j and makes it possible to extract the 



diffractive parton distribution functions of quarks and gluons |23l. |24||. Unlike the approach 
proposed by the Durham group, the Ingelman-Schlein (IS) proposal suggests the emission of 
Pomerons from the protons and, subsequently, the parton-parton interaction, as presented 
in Fig. [TJ As the Pomeron is a colorless object, the gluon emission leads to the Pomeron 
dissociation into hadronic states (X and Y in Fig. [1]), fragmenting into hadrons afterwards. 

Therefore, both processes have a similar structure but regard different phenomenological 
approaches to account for the diffractive interaction to produce the Higgs boson. In this 
work, we explore the diffractive factorization for the Higgs boson production at next-to- 
leading order (NLO), estimating the total cross section for the energy regime of the LHC. 
This study will allow us to confront our predictions with the ones obtained by the Durham 
group with the use of just a K factor to include all these corrections, since both approaches 
are evaluated at the same theoretical accurac}o So, we expect to see wether or not the curves 
1 It is important to note that the Durham group introduces the higher-order contributions to the gluon 
fusion vertex by a multiplicative factor of 1.5 [10(, as verified in Ref. 25]. It is not an accurate NLO 



calculation; however, we assume that the NLO contributions are included in that approach, allowing one 
to perform a direct comparison with our results. 



change their behavior with the inclusion of the NLO diagrams. Nevertheless, the addition 
of the GSP is mandatory in both cases in order to correctly predict the production cross 
section. As this probability is computed for a specific production mechanism in hadron- 
hadron collisions, the models for the survival factor employed here can be used in both 
mechanisms. Thus, this paper is organized as follows: in Sec. II, we present the inclusive 
production of the Higgs boson at NLO accuracy. Next, in Sec. Ill, the gluon-gluon luminosity 
is modified to describe the double Pomeron exchange between the protons. Further, in 
Sec. IV, the models for the GSP applied in this work are shown, as well as their estimations. 
In Sec. V, we discuss the main sources of uncertainties in this approach which affect our 
predictions. Then, in Sec. VI, our results are displayed as a function of the Higgs boson 
mass for different collider energies, as well as the rapidity distributions of the Higgs boson. 
Finally, in Sec. VII, we summarize our conclusions. 



II. INCLUSIVE PRODUCTION 



The inclusive process pp — > H + X at leading-order (LO) accuracy for the gluon fusion 
vertex corresponds to the fermionic triangular loop [2|, for which the quark top has the 
highest contribution 26]. Following the prescription of hard factorization given in Ref. 25], 
one is able to write the production cross section at LO as 

dC" 



plo{pp -> X + H + Y) = a r H 



dr 



H 



(1) 



where r# = Mfj / s is the Drell-Yan variable of the process, with s the center-of-mass energy 
squared, and do is a function of the variable tq = Mfj/Amg, defined as 



0o 



G F at 2 s (ii 2 ) 
288tt v / 2 



(2) 



with Gf = 1.16639 x 10~ 5 the Fermi constant, /i the renormalization scale, and Aq = 2[tq + 
(tq — 1)/(tq)]/tq. Considering the leading contribution of the top quark for the fermionic 
oop in the gg — >■ H vertex, we take the limit tq < 1 that corresponds to /(tq) = arcsin 2 Jtq 



25] . The last term in Eq.((T]) is the gluon-gluon luminosity 

dC" f 1 dx 



dr 



r 1 dr 

/ — g(x,M 2 )g(r/x,M 2 

Jr X 



(3) 



where M is the factorization scale, and g(x, M 2 ) is the integrated gluon distribution function. 
This luminosity will play an important role in this work, since it will be used to compute 
the NLO QCD corrections and will be modified to introduce the Pomeron structure function 
(PSF) to consider the emission of partons off the Pomeron in the double Pomeron exchange 
process. 

The singular and nonsingular virtual QCD corrections to the qq — > H vertex are expressed 
through the processes gg — > H(g), gq — > Hq, and qq — > Hg [25|, |27J]; the NLO contributions 
to the production cross section in pp collisions can be computed as follows: 



7T 



th-. h &v gg + Aa gq + Aa qq -, (4) 

dr H 



the singular virtual corrections being expressed in the function C(tq), and the nonsingular 
ones in the terms A<Tjj, with the renormalization and factorization scales fixed for the strong 
coupling constant a s (/i 2 ) and gluon distribution functions g(x,M 2 ), respectively. In order 
to add these corrections, the running strong coupling constant has to be input at NLO 
accuracy. We take into account the exact numerical solution of 



da s (/i 2 
din fi 2 



} 



167T 



(5) 



with fa = (11N C - 2N F )/3, 1 = {102N C - 38N F )/3, N c = 3, and N F depending on the 
quark flavor during the numerical calculation. 

The sin gula r virtual corrections correspond to the two-loop corrections, which are ex- 
pressed by [25] 



C(tq) 



7l~ 



+ <r Q ) + 
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2N F \ v 
log 



M 2 H 



(6) 



where n 2 refers to the infrared part of the cross section for real gluon emissions, and c(tq) = 
11/2 is solved analytically for tq = Mfj/Amq 27, 29)]. Next, the nonsingular virtual cor- 
rections Ao"y are obtained from the diagrams of gluon radiation in the gg and gq scattering, 
and qq annihilation. Each of them will be computed through Eq.(j3]), modified to include 



the q and q contributions 25] 
- 1 , &C™ a s (fi 2 



Act. 
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lofffl-f^ 
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-f P ss (f Jlog— + d gg {T, Tq) 



f[2-f(l-f)]log(l-f) 
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dr 



dC 9q a s (fx 2 



^ dr 
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1, M<1 i / 
-log— + log(l - r) 



^dC™ oJjj 2 ) , 



(7a) 

(7b) 
(7c) 



as f = th/t, and P gg {j) and P gq (f) are the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi split- 
ting functions [30J 
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1 -r. 
41 + (l-f) 2 



-r - 2 + r(l - r) V + 5(1 - r), 



(8a) 
(8b) 



3 f 

The F + is the plus distribution defined as F + (f) = F(f) — 5(1 — f) J Q df'F{f'). It is 
important to keep in mind the dependence of these individual cross sections on the parton- 
parton luminosities dC 11 / dr, since each of these luminosities will be replaced to introduce 
the diffractive interaction where the partons are be emitted off the Pomeron. 

As we are assuming the to p q uark contribution in this work, the dij(f, tq) functions can 



be evaluated numerically 27 



29|: 
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dgg{r,T Q ) = -y(l -f) 3 , 



d gq(T,T Q ) 



•l + 2f 



32 



(9a) 
(9b) 
(9c) 



Finally, we also account for the electroweak (EW) two-loop corrections to the gg — >■ H 



vertex 



31] , which increase the production cross section by 5% in comparison to the NNLO 



QCD corrections. Then, gathering all corrections to the inclusive Higgs boson production, 
the total cross section is written as 



Ctot — o"nlo(1 + <5ew)- 



(10) 
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Consequently, our predictions reproduce the results obtained in Ref. [25| for y/s = 14 TeV. 

However, even with a higher production cross section, the signal from the inclusive pro- 
duction is expected to be strongly attenuated by background processes. So, this is an op- 
portunity to explore the diffractive production as an alternative to detect the Higgs boson 
at the LHC. 



III. CENTRAL EXCLUSIVE DIFFRACTIVE PRODUCTION 



Our calculation is based on the IS approach, considering the Pomeron emission off the 
proton, followed by the interaction of its parton content that produces the Higgs boson, 
as illustrated by Fig. [TJ This procedure allows one to rewrite Eq.Q to replace the gluon 
densities by the Pomeron flux factor and its partonic distribution function. In this work, we 
assume a standard Pomeron flux, constrained from the experimental analysis of the diffrac- 



tive structure function performed by the HI Collaboration [24|. The PSF has been modeled 
in terms of a light flavor singlet distribution i.e., the u, d, and s quarks and their 

respective antiquarks. Also, a gluon distribution g(z) is included, with z the longitudinal 
momentum fraction of the parton in the hard subprocess. The Pomeron trajectory is as- 



sumed to be linear, a w (t) = ct^Cf 
to the data from the HI detector 



+ a'jpt, with a'jp and their uncertainties obtained from fits 



. We choose x^p where 
t m in\ ~ mpX%/(l — x w ) is the minimum kinematically accessible value of \t\, m p is the proton 
mass, and \t cut \ = 1 GeV 2 is the limit of the measurement. The HI parametrization provides 
two different inputs for the fit of the partonic structure functions. As our curves show very 
close results using both fits, we chose the fit A to perform our predictions, applying the cut 
x < x w < 0.05 in accord with this parametrization. 

Considering a Pomeron being emitted from each proton, one can express the luminosity 
of the CED process as 



r\r ij r 1 Ht Z" - 05 Hr 1 / t \ Z" 005 H-r- 

LU -CED _ aX UX IP ip I 1 x ,/2 / UX IP 77 f 2 




T 



—Fj/w/p Zp, -5—, M I ,(11) 



where ijr is the momentum fraction of the Pomeron carried by the parton relative to the 
hadron 1 (2), and the PSF 

Fi/w/p is expressed by 

Fi/w/ P = fw/p(xw)fi/ip ( — , M 2 ) , (12) 



fw/pi^Tp) being the Pomeron flux, and fi/-jp(/3, /i 2 ) the parton distribution function into the 
Pomeron, where i,j stands for g, q, and q. Nevertheless, it is necessary to correctly estimate 
the fraction of events during which the rapidity gaps will not be filled by the particles from 
underlying events with the use of the survival factor. 

In order to include the QCD and EW corrections to the gluon fusion vertex, as obtained 



in Ref. 



25 



31] ; we replace the parton luminosities of all individual cross sections, which 



keep the use of the MS factorization scheme. Moreover, we work in the same accuracy 
as the HI 2006 parametrization, which evolves the diffractive parton distribution functions 
with the DGLAP evolution equation at NLO accuracy |24j. The Pomeron-to-parton splitting 
functions, like Eqs. (l8"a|) and ( )8b|) . are not included in this evolution equation at this accuracy 
level. 

IV. GAP SURVIVAL PROBABILITY 



The computation of the total cross section for the CED Higgs boson production demands 
the addition of a survival factor, assumed enhanced diagrams for multi-Pomeron interactions 
in high energies [32] . The initial studies in this subject started computing the soft Pomeron 
exchanges through eikonal scattering and later added the enhanced diagrams with the triple 
Pomeron coupling which brings important contributions to the GSP. We adopt the 

two most discussed models for the GSP in CED Higgs boson production. 

Firstly, there is the Kaidalov-Khoze-Martin-Ryskin (KKMR) model that assumes 
enhanced diagrams for Pomeron exchanges, predicting a GSP for the CED Higgs boson 



production of 2.4% (1.5%) at 7 (14) TeV 



18 



. Secondly, there is the Gotsman-Levin-Maor 



(GLM) model 17|, working with the QCD and N=4 super Yang-Mills 19j, |33j. As a result, 



a probability was found o 
the GSP obtained in Refs. 



6, 



19 



(4.6%) at 7 (14) TeV. However, there is a change between 



33] . since the estimated GSP is increased from 0.15% to about 



3%. This effect is explained by the inclusion of semienhanced diagrams, that were neglected 
in previous calculations of the survival factor, and introduces important contributions. 

Moreover, there are other approaches that account for the GSP in CED Higgs boson 
production, indicating similar values as found in the KKMR model 34J but also very small 



2 Note that the probability of 1.5% in the KKMR model is seen as a lower limit for the GSP in the CED 
Higgs boson production at the LHC. This probability becomes higher depending on the Pomeron dynamics 



under consideration. 



s 



probabilities, like 0.44% 35j. Therefore, we employed the GSP from the KKMR and GLM 
models, currently in agreement with each other. 



V. RESULTS AND COMMENTS 



Evaluating Eq.(j4]) with the modified luminosity and including electroweak corrections, we 
have computed the cross section of the CED Higgs boson production for the LHC energies 
of 7 and 14 TeV. In order to include the suppression of the cross section due to underlying 
events, we have used the GSP for both kinematical regimes. In Figs. [2] and [31 we show 
our curves for the total cross section. Moreover, the cross sections using the KKMR and 
GLM models for the survival factor are also presented. In Fig. H] is displayed the rapidity 
distribution for a Higgs boson of M# = 120 GeV. As expected, the results lie in the same 
order as the ones from the Durham group. However, in comparison to our results for different 
LHC energies, we see a change in the shape of the curves, showing that the NLO diagrams, 
fully included in this approach, bring distinct contributions to the cross section. This aspect 
can not be predicted if one uses just a A-factor to include such corrections, which introduces 
further uncertainties to the predictions. 

Furthermore, the GSP is the only untested parameter for the exclusive Higgs boson 
production, causing an uncertainty of more than a factor of 2 if one compares the two GSPs 
applied to our predictions. As already successfully implemented for the Xc production in 
the Tevatron [7| , the estimations for the survival factor by both models seems to be reliable 
to be used in the CED Higgs boson production in the LHC kinematical regime. 

However, these small cross sections reveal that the detection of the Higgs boson in exclu- 
sive processes is going to be a challenge in the LHC, especially because the decay channel 
H — >■ bb 0]' with a higher branching ratio, is not going to be observed in the ATLAS 
and CMS experiments {38]. Then, other decay channels can be an alternative for this ob- 
servation, like H — > r + r~, but showing a lower branching ratio if compared to BR(66). In 
addition, maybe more feasible, one may look for the Higgs boson decay into a W + W~ pair 
in the range Mjj > 135 GeV, which has a higher branching ratio but smaller cross section 
than in the range Mh < 135 GeV, as seen in Figs. [2] and [3j 



VI. UNCERTAINTIES 



In the use of diffractive factorization, we found different sources of uncertainties that 
affect our results. First, the parametrization for the PSF provided by the HI Collaboration 
is constrained with the data obtained in HERA and does not contribute with significant 
uncertainties to the parton-parton luminosities. Next, the survival factor brings a large 
uncertainty to our calculations, since it is a model-dependent approach to account for the soft 
interactions between the protons. As presented in Sec. IIVt we apply two distinct approaches 
of the GSP for the production cross section, where the probabilities differ substantially 
between each other. Next, the perturbative calculations of the NLO corrections for the gluon 
fusion vertex have uncertainties related to the renormalization scale in the individual 
cross sections and also to the factorization scale fip, which mainly affects the parton-parton 
luminosities. 

In order to estimate the total effect of these uncertainties in our predictions, Fig. [5] 
shows the variation of the scales and the GSP with our results for y/s = 7 TeV. The upper 
curves are the predictions using the GLM model for the survival factor, with the dashed line 
representing the probabilities of 6.4%, with fi R = fip = M H . The dotted band shows the 
variation of the production cross section with the scales fiR and //p, where the lower limit 
represents the results with fin = fip = \Mh- Moreover, we performed the same analysis 
for the production cross section using the survival factor estimated by the Durham group. 
The results with a GSP of 2.4% are shown in Fig. [5] by the solid line, and the striped band 
represents the variation of the scales with the same limits as in the GLM results. 



As a result, the effect of all these parameters is an uncertainty of 5, which is in agreement 
with the uncertainties found in other approaches for the CED Higgs boson production 



The same analysis with sfs = 14 TeV leads to a similar effect on the predictions, showing an 
uncertainty close to 5. The gap seen in Fig. [5] between the results for the GLM and KKMR 
models is a result of our choice of the GSP from the KKMR model, since the probability of 
2.4% is taken as the lower bound for the survival factor in this approach for y/s = 7 TeV. 
In the case of higher values for the survival factor in the KKMR model, it is expected to see 
no gap between these curves. 
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VII. CONCLUSIONS 



In summary, we have computed the CED Higgs boson production at NLO accuracy 
for the kinematical regime of the LHC. In this calculation, we have included virtual QCD 
and electroweak corrections to the CED production of the Higgs boson in the diffractive 
factorization. As the nature of the underlying events is still an open question in high-energy 
physics, we have computed our predictions using two distinct models for the GSP in order to 
explore the possibilities of this production mechanism. Thus, we obtained a production cross 
section of 0.29 (1.2) fb at ^/i = 7 (14) TeV for M H = 120 GeV using the KKMR model and 
0.79 (3.68) fb at \fs = 7 (14) TeV in the GLM model. Our results agree with the ones from 
the Durham group, obtaining a production cross section of the order of a few femtobarns. 
Nevertheless, if the GSP is really of the order of 1%, the two-photon process becomes a 
oromising production mechanism in the LHC with a production cross section of 0.10-0.18 fb 
39[, with almost no suppression due to underlying events. Also, the cross section obtained 
in the photoproduction mechanism lies in the same range J40I, revealing that the scenario for 
the CED Higgs boson production is very competitive in the LHC energy regime. Therefore, 
both approaches for the CED Higgs boson production show similar results and may be an 
opportunity to discover the Higgs boson at the LHC, despite the fact that the observed 
decay channels will bring experimental difficulties for this detection. 
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FIG. 1. Feynman diagram representing the CED Higgs boson production in the IS approach for 
proton-proton collisions. The smaller blobs describe the parton content of the Pomeron, which 
results in the fusion of two gluons and the hadronic states X and Y. 
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FIG. 2. Total cross section in fb in function of the Higgs boson mass for the CED production at 
y/s = 7 TeV with no survival factor (solid line) and taking the GLM model (double-dot-dashed 
line) and the KKMR model (dashed line) for the GSP. 



16 




FIG. 3. The same as Fig. [2] for ^/i = 14 TeV. 
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FIG. 4. Differential cross section in fb in function of the Higgs boson rapidity for Mh = 120 GeV 
in different collider energies. The GSP is not included in this results. 
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FIG. 5. Production cross section in fb in function of the Higgs boson mass considering the vari- 
ation of the renormalization and factorization scales. The lines represent the predictions for the 
production cross section with hr = ptp = Mh using the survival factors from the GLM model (the 
dashed line for 3% and the dot-dashed line for 5%) and the one from the KKMR model (the solid 
line for 1.5%). Both dotted and striped bands represent the results with fin = fip = \Mu in the 
lower limit and fj,R = (j,p = AMh in the upper limit. 
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